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Abstract 

We perform certain alternating binomial summations with param- 
eters that occur in the analysis of algorithms. A combination of 
integral and special function and special number representations is 
used. The results are sufficiently general to subsume several previ- 
ously known cases. Extensions of the method are apparent and are 
outlined. 
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Introduction 



Alternating binomial sums arise frequently in computer science and data 
processing in the design and analysis of algorithms (e.g., P 1^ [TUl [TTl [T^ ) . 
The asymptotic form of such sums is often of interest in connection with 
determining the average- or worst-case run time. Because of the sign alter- 
nation of the summands, there may be substantial cancellation, masking the 
dominant behaviour. Recently there has been additional interest in certain 
alternating binomial sums and their connection with harmonic numbers and 
representation in terms of the Bell polynomials In fact, the ques- 

tion of a broader range of validity of such identities has been posed . In 
this article, we demonstrate a method to reach alternating binomial sum 
representations with a domain extended to the complex plane. 

Our approach is to obtain exact analytic relations. To these, known 
asymptotic relations may be applied if desired. Our results help to elu- 
cidate the connections between certain Bell polynomial representations, 
generalized harmonic numbers, Stirling numbers, and special values of the 
polygamma functions, and should be helpful in the analysis of either deter- 
ministic or probabilistic algorithms. Various integral representations pro- 
vide a convenient centerpoint of our development, but this is by no means 
necessary. There are many complementary approaches, including the use 
of finite difference operators [HEl. Afterall, A"/(.t) = {E - lYf{x) = 
(-l)"ELo(-l)'(fc)/(^ + ^here A/(a:) = f{x + 1) - f{x), Ef{x) = 
f{x + 1), and / is the identity operator. A way to think of the underlying 
combinatorics is in terms of number partitioning needed in the course of 
differentiating composite functions |19j . 

In this paper, we calculate the alternating binomial sums 

s(.,A.,™)^i:(^)(i^. (1) 

in multiple fashion for positive integers N and m and complex x £ C/{0, —1, 
—2, . . . , —N}. The special cases x = ±K for K a positive integer recover 
results of Kirschenhofer 8 and Larcombe et al. Jl]. The special case of 
a: = 1 in Eq. (1) has further applications in quantum information science 
[3 and this case is evident in our scheme. We point out additional special 



cases for x a rational number. 

Our work is more illustrative than exhaustive as there are many possible 
extensions. Especially when the summand in question contains a function 
for which Lemma 1 below applies, there will always be an equivalent rep- 
resentation in terms of the Bell polynomials Ym{xi, . . . ,Xm)- Moreover, 
as Eq. (10) below demonstrates, any time an expression contains Stirling 
numbers of the first kind, these may be replaced with Bell polynomials 
with generalized harmonic number arguments. Of the Bell polynomials, we 
note that they may be written as a lower triangular determinant save for a 
superdiagonal of — I's: 
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They satisfy the recursion relation 
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Reference JHl contains a background section on the Bell polynomials, wherein 
Proposition 2 repeats the well known determinant expression (2). Those 
authors also denote by "inverse Bell polynomials" that have usually 
been called logarithmic polynomials L„ (p. 140). Proposition 1 of Ref. 
|18) covers the parity of Bell polynomials when the even- or odd-indexed 
variables are put to zero and may instead by obtained from the determi- 
nantal expression (2). For further information on y„ we refer to standard 

works umini. 



Analytic number theory is an additional field where significant alternat- 
ing sums occur and these applications should not be overlooked. In fact, 
we take up an important such instance elsewhere. 



Summary of results and preparation 

Let TV > and m > be integers, x G C/{0, -1, -2, . . . , -iV}, V the 
Gamma function, (a)„ = V(a + n)/r(a) the Pochhammer symbol, B the 
Beta function, and ^Fq the generalized hypergeometric function [21 113| . 
Let s(n, fc) and fc) be Stirling numbers of the first and second kind, 
respectively P 12 El UHl [T7j . We have 
Proposition 1 
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InEq. (II), 
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g^^\x) = ^'^^'^(x) - i:^^^\x + + 1) = -(-1)^^! ^ 



fe=0 
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(126) 



where ijj = F'/F is the digamma function and ijA^"^ the polygamma func- 
tion The summation expressions in Eq. (12) follow by the use of the 
functional equations of these functions. In Eq. (10), Hn'^ are generalized 
harmonic numbers, 
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The proof of the equivalence given in Eq. (11) makes use of the following. 
Lemma 1. For differentiable functions / and g such that f'{x) = f{x)g{x), 
assuming all higher order derivatives exist, we have 



dx 



f{x)^f{x)Y, \g{x),g\x),...,g^^-^\x) 



(14) 



Proof of Lemma 1. Under the premise, g{x) = {d/ dx) \vl f{x) whenever 
/ ^ and {d/dxy f = /exp(— In f){d/dxy exp(ln/). The conclusion then 
follows as a special case of the Faa di Bruno formula for the derivative of 
a composite function. The result extends to x £ C when a branch cut for 
the logarithm is taken from the origin to the point at infinity. 

Remarks and an example, (i) The (exponential) complete Bell 
polynomials Yn may be obtained as a sum over the (exponential) partial 
BeU polynomials Bn,k- Yn{xi, . . . , Xn) = J2k=i ^nM^i^^2, ■ ■ ■ ,Xn-k+i), 
Yo — I, Yi{xi) — xi. (ii) The condition of the Lemma occurs often in 
practice and a nice example is given by the Gamma function with F' = Tip 
(m, p. 175). Then F(")(1) = r„(-7, a;2, X3, . . . , a;„) = J^" e"" In" x dx, 
where 7 = —"0(1) ^^^d xj = (—iy{j — 1)!C(.?) with ({s) the Riemann zeta 
function. 

(iii) The equality of Eqs. (3) and (4) is obvious. The original binomial 
series form (1) may be returned from Eq. (4) by the following two steps. 
First, depending upon whether N is even or odd we substitute into Eq. (4) 
either the expansion (^5_, p. 25) 
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or 

sinh^"-! X = y^(~l)"^^''2 sinh(2n - 2fc - l)a;. (156) 

k=0 ^ ^ 

We then apply tabulated integrals (0, P- 360) to find Eq. (1). 

(iv) The upper limit on the summation in Eq. (1) could just as well be 
put to cxD due to the property (^) — for k > n. 

(v) We are using notation for the Stirling numbers as followed by Comtet 
0] and Riordan and the reader should be aware of other conventions. 
Indeed the notation for these numbers has never been standardized ^ . (vi) 
The Stirling numbers of the second kind arc of rank one while of the first 
kind are of rank two. I.e., the latter numbers require a double summation 
in order to be expressed in terms of elementary factors 0j. (vii) We do 
not require them here, but mention that asymptotic forms of the Stirling 
numbers and functions are known. 

Proof of Proposition 1 

We now proceed systematically through the list of equivalences given in 
Proposition 1. In writing Eq. (1) in the form (4) we use the power series 
form of m+iFm [2] and apply {x)k/{x + l)k = x/{x + k). The terminating 
hypergeometric series in Eq. (4) is m + A^-balanced since the sums of 
numerator and denominator parameters differ by this positive integer. 

We next recognize that 

f^^\k)x + k r(a; + iV+l) {x)n+i x{x + 1)n ^' 

(16) 

that is equivalent to a partial fractional decomposition. Then 
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^{^l)"'-' t"'~'e-''\l-e-ydt. (17) 
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In the above the interchange of differentiation and integration is justified 
by the absolute convergence of the integral. We have shown the equality of 
Eq. (1) and (4)-(6). Equation (7) follows from (6) by using the definition 
of the hyperbolic sine function in terms of exponentials. 
To obtain the form of Eq. (8) we write Eq. (6) as 

S{x,N,m)^- — / f"~ie-("+^)*(e* - l)^dt, (18) 

(m - 1)! 7o 

and apply a generating function for the Stirling numbers of the second kind 

oo „ 

(e=^ = to! ^ S'(n,m)^. (19) 
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In obtaining Eq. (9) we first make the change of variable v{t) = 1 — e * in 
Eq. (6), giving 

C_-|\m-l /•! 



S(x,N,m) = ^^ — '—r v'" (I ~ vf-' In'^-'il - v) dv. (20) 
(m - 1)! Jo 

We then use a generating function for Stirling numbers of the first kind, 

°° x" 

In"^ (1 + x) = ml y s(n,m) — -, (21) 

n—m 

and carry out the integration with the Beta function. 

For Eq. (10) we apply Theorem B of Ch. V of Ref. U (p. 217) for the 
unsigned Stirling number of the first kind. We write this result in the form 

\sin + 1, A: + 1)1 = (-l)"+^s(n + 1, fc + 1) = yYk[H,,, ~Hi'\2\H^^\ 

{-l)'-\k-l)\Hi% (22) 

We then substitute for s{n,m — 1) in Eq. (9). 

In order to obtain the form (11) we use Eq. (5) and apply Lemma 1 
with the function f{x) — Nl/ {x)n+i, such that 

^f{x)^f{xmx)~i;ix + N + l)], (23) 



providing the function g(x) presented in Eq. (12a). Hence with this par- 
ticular g(x) we have 

(^)'/(:^)-/Wr, [g{x),g'{x),...,g^^''\x)] , (24) 

and the rest of Proposition 1 follows. 

Remarks, (i) There are many other variations on the possible generat- 
ing functions that may be introduced into the integrand of Eq. (6) to pro- 
duce equivalent forms with the Stirling numbers. Ref. |2] provides examples 
of these alternatives, (ii) The lower limit of summation in Eqs. (8) and (9) 
could just as well be put to due to the property s(n, k) = S(n, k) — 
for n < k. (iii) Equations (9) and (20) make it strikingly apparent how the 
case of a; = 1 is special, when the Beta function is no longer required, (iv) 
Single or repeated integration by parts in Eq. (6) permits the derivation 
of other forms of S{x, N, m) and of recursion relations for these sums. For 
instance, from e"^* = -{l/x){d/dt)e-'-'' and [e* /{N+l)]{d/dt){l-e-')'^+^ 
we immediately have 

S(x,N,m,) = -\S(x,N,m-l) + NS(x -1,N -l,m)], (25a) 

X 

and 

S{x,N,m) = j;j^[{x-l)S{x-l,N+l,m)~S{x~l,N+l,m-l)], (256) 

respectively. These relations have been obtained subject to Re a; > for 
Eq. (25a) and Re a; > 1 for Eq. (25b). 

Special cases 

We very briefly mention cases where x is an integer or a rational number 
in Eqs. (1) and (11). When a; = 1 we have for Eq. (12) 

N+l 

g(l) = V(l) - + 2) = - 5] - = -i/AT+i, (26a) 
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and 
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gi^l) - -(-1)^^! E - -i-lYm^^lW (266) 
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(r) 

where Hp are generalized harmonic numbers. These polygamma values 
are well known to relate to differences of the Riemann zeta function at 
integer argument, as V'^^Hl) = -(-1)'^^!C(^+1)- Here /(I) = iV!/(7V+l)! = 
l/(iV + 1) and then 



1 



x=l 

(27) 

When a; is a rational number it is possible to re-express the necessary 
derivatives g''^'\x). This can be done either in terms of the polygamma 
function or in terms of the Hurwitz and Riemann zeta functions (e.g., {12j). 
For instance we have ^(^)(l/2) = {-lY+^ £\{2^+^ - 1)C(^ + 1). 

When X is an integer we are able to express the sums of Eq. (12) in 
terms of generalized harmonic numbers. For instance we have for K a 
positive integer 



k = 



(28a) 
and 

g^'\K) = {-irH\[H''^Xi-Hp_'^ (285) 

Extensions 

The approach of this article may be extended to a great many other 
integrals. We outline some of this using the Beta function as the base. 
However, one could just as well apply the techniques to the confluent hy- 
pergeometric function i.F'i, the Gauss hypergeometric function 2-^1, and 
then to pFq more generally. 

We consider 

B{x,y)= u''-\l-u)y-^du = 2 sm'^'' ' ^ <j)cos^y- ^ cf) d(l) ^ B{y,x), 
Jo Jo 

min[Re x, Re y] > 0, (29) 

so that 

S{x,y,m,n) = [ u"" - \l - u)y-'^ In"'- ^ u\n''- \l ~ u) du 



Just as generalized binomial expansion gives 



4-Y (i) B{x,y). (30) 



B(.,,) = ^(-1F (^ /)-!-, (31) 



we have 

°° 1 / n \ 1 



5(a.,y,m,n) = (-l)'"-i(m-l)!X:-y(^j (l-y).^^^- (32) 



In obtaining this equation we wrote j^) ~ ("1)"'(1 ^ y)j/j- Of course 



3=0 ■ 

the sum terminates when y — 1 is a positive integer and we again exclude 
nonpositive integer values for x. The Pochhammer polynomial occurring in 
Eq. (32) has derivative 

^(1 - y)j = (1 - vUi^il -y)- + 1 - y)]. (33) 

Therefore, for instance, Lemma 1 may be used in finding all higher order 
derivatives. 

We may next change variable in Eq. (30) to obtain 

POO 

S{x,y,m,n) = {-ly-^ I w"-Mn™-i(l-e-'')(l-e-'')^-ie-^''rft;. (34) 
Jo 

Otherwise, we may first put 2; = 1/u in Eq. (29) and v = \nz, giving 

B{x,y)= / e-^'=+y-'^^''{e'" -ly-^dv. (35) 
Jo 

Then we have 

/•OO 

S{x,y,m,n) = / e-(^+2'-^>(e'' - 1)^"^ [(-l)'"-i?;'"-Mn"-i(e" - 1) 
Jo 

Equations (34) and (36) are forms suitable for re-expression as a summation 
with Stirling number coefficients. 

By the same token, one may repeatedly integrate the Beta function, 
thereby obtaining binomial summation expressions for integrals of the form 

/•I n _ uPYuy~-^ 

I{x,y,p,m)= du, (37) 

^0 In u 



for Re a; > 0, Re y > 0, and m an integer. In this way, we obtain extensions 
of tabulated integrals such as given in Sections 4.267 and 4.268 of Ref. 

Final remarks. We could also develop series representations with Stir- 
ling number coefficients using divided difference formulas. For example, we 
have P 



these formulas also exhibiting the inverse relations possible with the Stirling 
numbers. Especially for the Gamma and hence the Beta function, the finite 
differences and derivatives are relatively easily determined, due to their 
respective functional equations. 
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